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ABSTRACT 

We present a unifying scenario to address the physical origin of the diversity of X-ray lightcurves 
within the p variabihty class of the microquasar GRS 1915+105. This 'heartbeat' state is characterized 
by a bright flare that recurs every ^ 50 — 100 seconds, but the profile and duration of the flares varies 
significantly from observation to observation. Based on a comprehensive, phase-resolved study of 
heartbeats in the Rossi X-ray Timing Explorer archive, we demonstrate that very different X-ray 
lightcurves do not require origins in different accretion processes. Indeed, our detailed comparison 
of the phase-resolved spectra of a double-peaked oscillation and a single-peaked oscillation shows 
that different cycles can have basically similar X-ray spectral evolution. We argue that all heartbeat 
oscillations can be understood as the result of a combination of a thermal-viscous radiation pressure 
instability, a local Eddington limit in the disk, and a sudden, radiation-pressure-driven evaporation or 
ejection event in the inner accretion disk. This ejection appears to be a universal, fundamental part 
of the p state, and is largely responsible for a hard X-ray pulse seen in the lightcurve of all cycles. We 
suggest that the detailed shape of oscillations in the mass accretion rate through the disk is responsible 
for the phenomenological differences between different p-type lightcurves, and we discuss how future 
time-dependent simulations of disk instabilities may provide new insights into the role of radiation 
pressure in the accretion fiow. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: individual (GRS 
1915+105) — X-rays: binaries 



1. INTRODUCTION 

As one of the brightest and most variable sources in 
the X-ray sky, the microquasar GRS 1915+105 is an ex- 
cellent case study for the investigat ion of evolving accre- 
tion flows. Discovered in 1992 by ICastro-Tirado et al.l . 
the black hole has been in outburst for over 17 
years and has b ecome famous for its superluminal jets 
(iMirabel &: Rod riguez 1994), its bizar re X-ray variabil- 
itv (IBel loni et al...2000. hereaf ter BOO; IKlein-Wolt et al.l 
120021 : iHannika inen et all|2005[ ), and its disk-jet interac- 
tions (Fender & Belloni' '2004'; 'Fender. Belloni. fc Gallo 
2004; Ncilscn & Lee 2009; Neilsen. Rcmillard, fc Lee 
2dllaD 



Of its many classes of X-ray variability, the best stud- 
ied are the x st ate, which produces steady optically 
thick jets (see, e.g. lDhawan et al.]l2000l : iKlein-Wolt et al.l 
l2002f ). the /? sta te, a wild 30-minute cycle with discrete 
ejection events (jMirabel et al.l [19981) . and the p state, 
which is affectionately known as the 'heartbeat' state for 
the similarity of its lightcurve to an electrocardiogram 
(see Figure [T]) . The p-type cycle consists of a slow rise 
followed by a sh ort bright pulse, repeating with a p eriod 
of roughly 50 s (|Taam. Chen, fc Swanklll997l :lB?)m. 

After nearly two decades of X-ray monitoring of 
GRS 1915+105, it is clear that its remarkable diver- 
sity of X-ray states is also reflected in the states them- 
selves. Indeed, ev en in their original paper on these 
variability classes, iBOd divide the x class into four 
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Fig. 1. — Sample cycles from PCA observations of single-peaked 
(left) and double-peaked (right) oscillations: (top) count rate 
lightcurves, (middle) HRl, and (bottom) HR2. These sample cy- 
cles are taken from observations 40703-01-07-00 and 60405-01-02- 
00, respectively. In this paper we analyze all RXTE observations of 
p-like cycles; the data here are representative of the broad behavior. 
For reference, we have marked peak count rate times (identified by 
cross-correlation; see Section |2]l with dashed vertical lines. 



sub-categories that differ in hardness and noise prop- 
erties. Furthermore, some observations may exhibit 
characteristics of multiple states (e.g. /3-type variations; 
iNeilsen. Petschek. fc Leell2011b[ ). and various instances 
of a single state may even differ in their l arge-sc al e vari- 
ability patterns. For example, Massar o et al.l (J20101 
hereafter MIO) recently analyzed a long iJeppo SAX ob- 
servation of GRS 1915+105 in the p state and reported 
significant changes in the lightcurve shape: the number of 
pulses or peaks per cycle ranged from one to at least four. 
These pulses can be separated by 10 seconds or more, and 
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Fig. 2. — (Left) The measured period of the p cycle for all the analyzed observations. Periods cluster around 50-60 seconds, although 
there are some excursions as high as 340 seconds. The inset shows the drift and scatter that characterize the cycle period in a typical 
observation. (Right) The distribution of periods colored according to peak multiplicity. 



the later peaks are typically harder (jTaam et al.l 119971 : 
IPaul et all [19981) . 

In the hopes of uncovering direct evidence for the 
origin of these luminous p-type pulses, we performed 
the first joint Chandra/RXTE phase-resolved spectral 
analysis of an observation of GRS 1915-1-105 in the p 
state (Ncilscn ct al. 2011a). Using the Chandra grat- 
ings (jCanizares et al.l |2005j), we showed for the first 
time that changes in the broadband X-ray spectrum 
drive physical changes in the accretion disk wind on 
time scales as short as 5 seconds, and that this wind 
is sufficiently massive to cause state transitions in the 
disk. Based on RXTE spectra, we argued that ra- 
diation pressure plays several key roles in the accre- 
tion disk, from driving the observed limit cycle via 
the Lightman-Eardley instability (a.k.a. radiation pres- 
sure i nstability: Lightm an fc Eardlevlll974l : iBelloni et al.l 
Il997t iJaniuk et al.i |20Q0( ) to literally pushing the inner 
edge of the accretion disk away from the black hole (a 
local Eddington limit: iF ukuc 2004; Heinzeller fc DuschJI 
I2007t iLin. Remillard. fc Homan .20091) until the global 
Eddington limit is reached and the disk is evaporated 
or ejected at the ma ximum accretion rate (see also 
IJaniuk fc Czernvl 120051 ). For this last point, we argued 
that the second peak in our double-peaked lightcurves 
could be explained as bremsstrahlung emitted when the 
ejected gas collides with the corona. 

But what is the significance of these results whe n som e 
instances of the p state only have a single peak ([MIOT )? 
That is, if there is no "second peak" in the lightcurve, 
is mass ejection still required to explain the X-ray os- 
cillation? Are we to understand different variations of 
the p-like cycles as driven by the same physical processes 
or as producing similar lightcurves by coincidence? In 
this paper, we begin to address the remarkable diversity 
within the variability classes of GRS 1915-1-105, with a 
comprehensive phase-resolved study of all the heartbeat- 
like states in the RXTE archive. 

We describe the observations, data analysis, and phase- 



folding in Section 2. In Section 3, we perform a detailed 
comparison of the timing and spectral behavior of single- 
and double-peaked p states. In Section 4, we interpret 
our results as clear evidence that a single physical sce- 
nario may produce a wide variety of lightcurves, and we 
consider the e xact origin of the variations observed by 
iMassaro et all (2010). 

2. OBSERVATIONS AND PHASE FOLDING 



Since its launch in 1995 (|Jahoda et al.lll996[ l. RXTE 
has observed GRS 1915+105 roughly twice a week, for a 
total of over 1600 pointings as of the original writing of 
this paper. In this paper, we will focus exclusively on the 
data from the Proportional Counter Array (PCA), which 
covers approximately 2-60 keV. We extract 1-second nor- 
malized lightcurves from the binned data modes in three 
energy channels: A = 2 — 5 keV, B = 5.2 — 12 keV, and 
C = 12 — 45 keV. The lightcurves are normalized using 
the PCA lightcurve of the Crab nebula. That is, for each 
PCU, we renormalize the mission-long Crab lightcurve to 
count rate values 1100, 1140, and 330 counts per second 
for channels A, B, and C, respectively. This roughly 
corresponds to the raw count rate of PCU2 around MJD 
52000. We then apply these same normalizations to the 
GRS 1915-1-105 count rates so that the intensity scale 
is the same for all observations. We also produce two 
hardness ratios: HR1= B/A and HR2=C/S. 

In order to identify p state observations, we examined 
1-second lightcurves and color-color diagrams (CD) for 
all observations of GRS 1915-1-105. We selected observa- 
tions with regular or quasi- regular bursts of the appropri- 
ate shape, duration, and CD by visual inspection. This 
task can be rather tricky since there are variations within 
the heartbeat state (Section [T]) as well as strong resem- 
blances between the p class and some portions of the v 
and K classes; our analysis may include some cycles from 
these classes. However, the fact that our phase-folding 
method can be used for all these quasi-regular cycles is 
indicative of the physical robustness of our results. We 
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Fig. 3. — Average phase- folded lightcurves and hardness ratios for the single- and double-peaked oscillations depicted in Figure [T] i.e. 
RXTE observations observations 40703-01-07-00 and 60405-01-02-00, respectively. Two cycles are shown for clarity. The smooth curves 
represent the average lightcurve profiles, and the data are plotted in gray (204 cycles on the left, 265 cycles on the right). It is obvious that 
despite the scatter and variations in the period, the phase profile of the p oscillation is remarkably constant from cycle to cycle. Note that 
the energy ranges plotted here differ from iNeilsen et al.i (i2Clllai ) because we are using new, normalized lightcurves (Section [2]l. 



show samples of the Ughtcurves and hardness ratios for 
two typical (long) observations in Figure [T] {RXTE ob- 
servations 40703-01-02-00 and 60405-01-02-00 on the left 
and right, respectively). The oscillation is obvious in 
all panels. Inspection of these data and our ensemble 
of lightcurves reveals the same variations noted by MIO, 
namely that the number of strong peaks in the lightcurve 
changes over time scales of days or longer, as do the du- 
rations of the pulses and the time delays between them. 
As we will see later that there is likely a continuum of 
lightcurve shapes, we will not devote any time here to 
detailed phenomenological categorizations of these pro- 
files.^ 

In iNeilsen et al.l ()2011al ) , we showed that we can very 
accurately characterize the physics of this oscillation by 
phase-folding the individual cycles, which takes into ac- 
count the variable oscillation period. In other words, we 
stretch or compress and then combine all the individ- 
ual cycles in a given observation, and study the spectral 
conditions at each part of the cycle. We determine the 
start times of each cycle by means of an iterative cross 
correlation method t hat is described in more detail in 
INeilsen et al.l ()2011al ). Briefly, we use cross-correlations 
to identify the main peak in the count rate for each cy- 
cle, and then define the times of peak count rate as phase 
= 0. With this analysis, we measure 10068 peak times 



in 242 observations. For reference, we plot the resulting 
cycle periods (i.e. the distance between successive times 
of = 0) in Figure [2j It is evident from this figure that 
typical periods are < 100 seconds, although some single 
cycles may last over 300 seconds. Within a given obser- 
vation, there is a fair amount (> 10%) of cycle-to-cycle 
scatter in the period, but the period is usually very stable 
around the mean (see inset. Fig. [5]). 

Once we have defined a phase ephemeris for all 242 ob- 
servations, we create phase-folded lightcurves for each, 
and then we identify them as single-peaked, double- 
peaked, or as intermediate cases. Observations may be 
classified as intermediate if there is substantial noise (e.g. 
if the observation only includes a few cycles) or if the 
cycles have unusual shapes or strong period variability 
(like K-type cycles), so that it is difficult to reliably dis- 
tinguish two individual pulses from a single pulse. Of the 
242 observations included here, we classify 101 as hav- 
ing double-peaked cycles, 77 as single-peaked, and the 
remaining 64 as intermediate cases. Because there is an 
inherent difficulty in deciding whether structure in the 
lightcurve is noise or a physical peak, there is an uncer- 
tainty in these numbers of about 10-15 (estimated by 
repeating the classification). This uncertainty has a neg- 
ligible effect on the numbers presented in this section and 
no effect on the subsequent spectroscopy our physical in- 
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terpretation. 

There are some statistical differences in the period be- 
tween cycles with different numbers of peaks. The right 
panel of Figure [2] shows that single- and double-peaked 
cycles have similar periods, at 61 ± 23 s and 64 ± 18 
s, respectively, while intermediate cycles are generally 
shorter, with an overall mean of 52 ± 27 s. The uncer- 
tainties represent the sample standard deviations; the 
standard errors of the mean are < 0.5 s. At very short 
periods, the intermediate cycles show a bimodal period 
distribution, which is reminiscent of the k stat e (BOO ) 
or the irregular bursts studied by lYadav et al.l (|1999l ). 
However, we will see that these intermediate cases exhibit 
similar spectral behavior to single- and double-peaked cy- 
cle types, so they are unlikely to represent a completely 
different physical phenomenon. 

Figure |3] shows the phase-folded lightcurves for the 
single- and double-peaked observations illustrated in Fig- 
ure [TJ The average phase-folded lightcurve profiles, 
shown as solid lines, are superimposed on the individ- 
ual cycles from these same observations shown as gray 
points. There are a total of 204 cycles and 265 cycles 
on the left and right, respectively. The scatter is limited 
and it is clear that within an observation the individual 
cycles have the same profile. On the other hand, it is 
also clear that heartbeats can differ substantially from 
one observation to the next (i.e. in terms of the num- 
ber of peaks or pulses ). Building on our p rior study of 
double-peaked cycles (Neil sen et al.l[2011al ). we focus in 
this paper on a comparison between single-peaked cycles 
and double-peaked cycles, and the question: can these 
two cycle types be produced by the same mechanism? 
In other words, does the single pulse correspond to the 
first (soft) pulse of a double-peaked cycle, the second 
(hard) pulse, a combination of the two pulses, or is it 
completely different? 

3. THE SINGLE PULSE 

In this section, based on timing analysis and broad- 
band spectroscopy, we develop a case that the single pulse 
in single-peaked cycles corresponds to the second, hard 
pulse in double-peaked cycles. As such, our results will 
indicate that the physical processes responsible for the 
production of this hard pulse are active in all p-like cy- 
cles. In other words, these processes are fundamental to 
the heartbeat state. 

3.1. Folded Lightcurves and Timing Analysis 

The first indications that the single pulse is analogous 
to the hard pulse come from the lightcurves and phase- 
folded lightcurves in Figures [1] and 111 In these plots, HR2 
rises sharply from its minimum value around = 0, giv- 
ing the appearance of a local maximum. The local max- 
imum in HR2 coincides with the maximum value of HRl 
in both single-peaked and double-peaked cycles. This 
phase of "maximum spectral hardness" occurs during the 
main pulse in single-peaked cycles, while it is very close 
to the second pulse when there are two peaks in the cycle. 

An alternative indication of this same phenomenon is 
the delay between various maxima in the hardness ratios 
and in the hard and soft X-ray lightcurves. For example, 
consider that in a double-peaked cycle, HRl peaks at 
the same time as the C-band lightcurve, but after the 
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Fig. 4. — Distribution of phase delays between maxima in (a) the 
A and C bands and (b) the C band and HRl, for all observations 
analyzed in this paper. The top panel shows that single-peaked 
cycles generally have a small delay between their soft and hard 
maxima, but double-peaked oscillations have a much longer delay. 
In this respect, the intermediate cases sometimes mimic single- 
pulse profiles and sometimes behave like double-peaked cycles. The 
bottom panel demonstrates that there is no significant difference 
between any cycle types in the delay between the maximum in 
HRl and the maximum in the hard X-ray lightcurve, although HRl 
typically peaks a little later than the C-band lightcurve itself. 



A band, while all three quantities are synchronized in 
single-peaked observations. 

The delays between these maxima are easily quanti- 
fied from the phase-folded lightcurves, and we present 
the phase delay distributions in Figure SI It should come 
as no surprise that double-peaked oscillations exhibit a 
significant delay between the hard and soft X-ray pulses 
(top panel): A0c^^ = 0.08 ± 0.05, where the uncer- 
tainty is the sample standard deviation. As the average 
period of these oscillations is ^ 64 seconds, we find that 
the hard pulse typically follows the soft pulse by ~ 4.9 
seconds. For single-peaked cycles, the delay between the 
peaks in the A and C bands is much smaller (consistent 
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Fig. 5. — The distribution of HR2 at two phases of the heartbeat 
cycle for all observations analyzed in this paper: (top) (/> = 0, and 
(bottom) at the maximum of the C-band lightcurve. The top panel 
shows that double-peaked cycles are typically softer than single- 
peaked cycles at their respective times of 4> = 0, while the bottom 
panel shows that both types of cycle have very similar hardness 
near the maximum of the C-band lightcurve. Thus the single pulse 
may correspond to the hard (second) pulse of double-peaked cycles. 

with zero): A(j)c-A = 0.02±0.02. The intermediate cases 
show a bimodal distribution, apparently including mixed 
contributions from single-peaked and double-peaked cy- 
cles. As suggested above, all cycle types exhibit a very 
short delay between the maximum of HRl and the hard 
X-ray pulse. This delay is A^hri-c = 0.02 ± 0.04 for 
single-peaked cycles, 0.01 ±0.04 for double-peaked cycles, 
and 0.01 ± 0.03 for intermediate cases. 

Also relevant are the values of the hardness ratios 
themselves. In Figure [5] we present the distribution of 
HR2 at (j) — 0, when the count rate is at a maximum 
(top), and at the maximum of the C-band lightcurve 
(bottom). In general, we see that at = 0, double- 
peaked oscillations are softer than cycles with a single 
peak (intermediate cases lie in between the two). This 
is to be expected if the single pulse corresponds to the 
hard pulse of a double-peaked cycle. That is, when there 



are two peaks, our cross-correlation method picks out 
the brighter, softer first pulse for (j> = 0; in the absence 
of the soft pulse (i.e. one pulse only), cj) = coincides 
with the hard pulse. The bottom panel of Figure [5] con- 
firms this conjecture: all cycles have very similar HR2 
when the C-band flux is maximized. This suggests that 
a common light source produces a pulse of hard X-rays 
in every cycle. 

Again, the relative timing of the X-ray flux and spec- 
tral hardness ratios indicates that the intermediate cases 
represent a continuum of profiles between single- and 
double-peaked oscillations. A completely different physi- 
cal process might be less likely to produce the same tim- 
ing signatures, so we conclude that at 1-second resolu- 
tion, the variations in the cycle profile can be explained 
in the context of two main pulses, one soft and one hard, 
whose relative strengths, duration, and separation can 
vary from observation to observation. In Section 13.21 we 
employ phase-resolved spectra to unify these pulses in a 
single physical framework. 

3.2. Phase- Resolved Spectroscopy 

In lNeilsen et al.l (j2011af) . we showed that in a double- 
peaked cycle {RXTE ObsID 60405-01-02-00), the X-ray 
spectrum of the hard pulse was described very well as the 
sum of two major components: a disk blackbody and 
a power-law type component, possibly bremsstrahlung 
or Comptonized disk emission. Statistically, we could 
not distinguish between different Comptonization models 
and bremsstrahlung emission, but all models required a 
sudden decrease in the electron temperature associated 
with this component during the hard pulse. 

In this section, we ask whether an observation of a 
single-peaked cycle can be described by the same model 
with similar phase dependence; we return to the physi- 
cal interpretation of this model in Section U) We per- 
form phase-resolved spectroscopy of the single-peaked 
observation represented in Figures [T] and |31 RXTE Ob- 
sID 40703-01-07-00, and co mpare to our re sults on the 
double-peaked cycle in Neils en et al.l ()2011aD . We implic- 
itly assume that these two observations are representa- 
tive of single- and double-peaked heartbeats. This as- 
sumption is justified by our results in Section [XTl which 
demonstrate clear similarities between different observa- 
tions of a given type of cycle. 

Our spectral model consi sts of interstellar absorpt ion 
by cold gas ftbabs: Wilms. Allen, fc McCravl 120001 ) . a 
hot disk (ezdiskbb; [Zimmerman et a l. 2005) convolve d 
through a scattering kernel (simpl: [Steiner et alJl2009D . 
a Gaussian emission line at 6.4 keV, and a high-energy 
cutoff (highecut, with functional form exp(— .E/E'foid) 
for E > Ecut', we set E'cut — 0). simpl takes a seed spec- 
trum and scatters a fraction fsc of the source photons 
into a power law, approximating the high-temperature, 
low optical depth regime of Comptonization. We use 
highecut to account for curvature in the hard X-ray 
spectrum, and we prefer simpl to powerlaw because 
it conserves photons and includes a physically-realistic 
rollover at low energies. Our method is precisely the 
same as that reported in INeilsen et all (j2011aD : we 
extract phase-resolved spectra and responses and, ap- 
plying the model described above, we fit the emission 
from 3-45 keV. A ll spectral fitting is done in ISIS 
(jHouck fc Denicolal ,200ft Houck, ; 2002,) . We assume a 
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Fig. 6. — Model fit parameters and their 90% confidence limits as a function of phase for the PC A spectra of the single-peaked heartbeat 
oscillation (filled black points; RXTE observation 40703-01-07-00). For comparison, we overplot our parameters for a double-peaked cycle 
(open red circles; Noilscn ot al. 2011a, RXTE observation 60405-01-02-00). Two cycles are shown for clarity. Since the oscillation periods 
are near 50 seconds, each phase bin corresponds to roughly one second of real time. The fits are extremely similar; in conjunction with the 
phase shift between the two parameter sets, this result implies that the pulse in single-peaked cycles physically corresponds to the hard 
second pulse in double-peaked cycles. 

Given the differences in the Ughtcurves (Figure [1]) , the 
similarity of the spectral fits is striking. The most im- 
portant difference is a clear phase shift: a given feature 
in a given parameter occurs later in phase in double- 
peaked cycles. In terms of the parameters themselves, 
the double-peaked cycle has a slightly higher and flatter 
baseline in Tobs, and the spike is smaller; the disk ra- 
dius increases for a longer interval relative to the single- 
peaked cycle. The double-peaked cycle also has a slightly 
larger photon index, scattering fraction, and electron 
temperature than the single-peaked oscillation. 

Two points from our spectral fits deserve special con- 
sideration: 

1. It is remarkable that in these two very differ- 
ent cycles, the measured radius of the disk plum- 
mets sharply and comes to rest (albeit briefly) 
at the same value of Rm ^ 40 km. The pos- 
sibility of a "touchdown" radius is especially in- 
teresting given that 40 km is the radius where 
one expects the peak emission from a thin disk 
around a maximally-spinning 14 Mq black hole 
^McClintock ct al.,,2006, ). 



dista nce and incl i natio n of Z? = 11.2 kpc and i — 
66UFender et alj[l99l); we fix iVn = 5 x 10^^ cm'^ 
(|Lee et al. 200^ and references therein) . Overall, the 
model provides an excellent description of the data, with 
a reduced xV^^ = 1847.6/2149 = 0.86. 

The resulting fit parameter values for our single-peaked 
oscillation are shown in black in Figure |6l The maxi- 
mum observed temperature in the disk, Tobs, is mostly 
constant near 1.1 kcV throughout the cycle but spikes to 
~ 2.2 keV at (/) = 0. While the disk temperature hovers 
around 1.1 keV, the inner disk radius i?in rises steadily 
from ~ 60 km to '^ 95 km, but then turns over and plum- 
mets quickly to ^ 40 km (we use a color-correction factor 
/ = 1.9, appropriate for high luminosity; J. Steiner, pri- 
vate communication). The simpl photon index F, the 
scattering fraction /sc, and the high-energy cutoff ii^foid 
all decrease steadily for most of the cycle, from roughly 4 
to 1.5, 0.9 to 0.1, and > 500 keV to 30 keV, respectively. 
Near (/) = 0, F and /sc rise sharply, but -Efoid (i-e. the 
electron temperature) dips to 5 keV. During this brief 
dip, there is also weak evidence of a coincident dip in F . 

For comparison, we overplot our (jNeilsen et aLll2011aD 
measured parameters for the double-peaked oscillation. 
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2. Both single-peaked and double-peaked oscillations 
are marked by a short interval where the scattering 
fraction approaches 1 and ii^foid drops to '^ 5 keV, 
indicating the sudden appearance of an optically- 
thick population of relatively cool electrons. In 
double-peaked cycles, this interval corresponds to 
the hard pulse. In single-peaked cycles, it coincides 
with the main peak in the lightcurve at ^ = 0. This 
connection confirms the identification of the single 
pulse as the hard pulse. 

On the whole, it is abundantly clear that the same spec- 
tral processes are operating in both of these cycles. 

4. DISCUSSION 

In Section [Sj our phase-resolved timing and spec- 
tral analysis demonstrated that single-peaked heartbeats 
should be interpreted as double-peaked cycles with a 
weak soft pulse, while the brief pulse of hard X-rays is 
truly fundamental to the p cycle. In this section, we con- 
sider the origin of these pulses and the implications of 
their relative variations. 

In our detailed phase-reso lved spectral analys is of 
a double-peaked oscillation (jNeilsen et al.l l2011a[ ). we 
explored two Comptonization models (simpl and 
nthcomp; Zdz iarski etall 119961: iZvcki et al] [1999) and 
also bremsstrahlung as descriptors of the broadband X- 
ray emission in the heartbeat state. While these models 
differ in their physical interpretations, they implied very 
similar behavior in the inner accretion disk, and all three 
models required the sudden appearance of a new popu- 
lation of electrons with a temperature of ~ 5 keV during 
the hard pulse. We argued that this was probably due to 
the sudden ejec tion of matter from the disk at very high 
luminosity (see iJaniuk fc Czernvl 120051 for a theoretical 
basis). 

To produce the waves in M that cause the system to 
exhibit these changes on a cyclic basis, we invoked a 
globa l Eddington instability fLightman fc E ardlevI 
119741 : iNeilsen et all [2011a). " This instabil ity, a 
class of thermal-viscous instability (iTaam et al.l 
Il997t iNavakshin. Rappaport. fc Melial l2000f l. has 
also been c alled the radiation pressure in stabil- 
ity (RPI: IJaniuk. Czernv. fc Siemiginowskal 120001 : 
IJaniuk fc Czernvl 120051 ). a term we adopt in what 
follows. We also found a long interval of increasing 
disk radius at constant temperature, which we argued 
was an indication of a local Eddingtoii limit in the dis k 
(|Fukuell2004l: IHeinzeller fc Dusch]|l2007l : ILin et al.ll2009l) . 
Local Eddington effects can occur in geometrically thin 
disks when the local vertical component of radiation 
pressure becomes sufhciently strong to overcome gravity, 
disrupting the disk interior to some critical radius. In 
short, when the RPI raises the local accretion rate, 
the disk responds via the local Eddington process: it 
maintains a constant temperature but increa ses its 
inner radius as it expels a ring of gas at R\y, (jFukud 
I2OOI IHeinzeller fc Duschll 120071 : JLin et al.l 1200 91 These 
phenomena (mass ejection and the local Eddington 
process) reflect the critical importance of radiation 
pressure in the accretion dynamics of double-peaked 
heartbeat states. 

Furthermore, in the single-peaked cycle we have stud- 
ied here, the same two effects (i.e. mass ejection and disk 
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Fig. 7. — Accretion rates versus phase implied by our fits to the 
phase-resolved PCA spectra of single- (black) and double-peaked 
heartbeats (red) . For direct comparison, we also overplot M for the 
double-peaked cycle, normalized to the mean value for the single- 
peaked cycle (green dashed line). In double-peaked cycles, the 
accretion rate is higher and the wave of matter is sustained for 
longer; in single-peaked cycles, the overall accretion rate is lower 
and spikes sharply around (/> = 0. We discuss in Section |4] how the 
differences in the M profile may be responsible for the different 
X-ray lightcurves. 



evolution via a local Eddington limit) also dominate the 
spectral evolution. Thus we believe that a single unified 
model can successfully describe the various lightcurves 
of p-like states, and we identify radiation pressure as the 
driver of this unified model. From the oscillating mass 
accretion rate (the radiation pressure instability) to the 
slow growth of the disk radius (the response of a luminous 
thin disk to the rising accretion rate) to the ejection or 
evaporation of the inner disk (a consequence of the sud- 
den final influx of matter), radiation pressure mediates 
all major aspects of the p cycle. 

If a common set of accretion processes can describe the 
range of p-like cycles seen in the RXTE archive, we must 
also address the question: why are some oscillations dif- 
ferent than others? If the hard pulse is always observed 
in the heartbeat state, where does the soft pulse go? To 
answer this important question, we look to the compar- 
ison of our spectral parameters in Figure jB] The salient 
features of the phase evolution of the disk are the plunge 
of i?in towards the black hole and the rapid heating of 
the inner disk. Ultimately, it is the timing of this plunge 
relative to the spike in the disk temperature that deter- 
mines the cycle profile and the relative strengths of the 
two pulses. 

But examining the evolution of i?in and Tobs in light of 
the role of radiation pressure in the heartbeat state, we 
see two puzzles that merit further consideration. First, 
why are the changes in the disk so catastrophic? In prin- 
ciple, the local Eddington process is reversible, so that 
i?in could also decrease slowly at roughly constant tem- 
perature (e.g. reversible changes tied to the local Edding- 
ton limit appear to operate along the flaring branch of 
Z-sources: ,Lin et al.. ,2009). But in the heartbeat state of 
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GRS 1915+105, the symmetry appears to be broken: the 
disk expands and contracts along vastly different paths 
in the i?in — Tobs plane. It is unclear whether this asym- 
metry is related to mass ejection, a steep gradient in the 
accretion rate, reaching some energetic hmit in the disk, 
or some other unknown factors, but it is clear that there 
is significantly more to the behavior of the disk than the 
local Eddington process. 

This brings us to the second puzzle: what causes the 
local Eddi ngton effect to stop acting on the disk? In 
iNeilsen et al. (2011a), we found that in a double-peaked 
cycle, i?in grows at roughly constant temperature until 
4> ~ 0.9, when the catastrophic changes take over to be- 
gin a new cycle in the disk. But for the single-peaked 
cycle shown in Figure [6l Rin stops growing and actu- 
ally decreases slightly from (j> ~ 0.75 — 0.9. This mo- 
mentary plateau in Ri^ signals the end of the local Ed- 
dington process, but here it occurs almost ten seconds 
before the catastrophic changes in the disk. In other 
words, the plunge of Rin does not cause the end of the 
local Eddington effect (although they may coincide, as in 
double-peaked cvcles: INeilsen et al1l2011a[ ). Instead, the 
end may be controlled by changes in the accretion rate. 
From our fits, this puzzling phase interval shows a steady 
Rin plus an increasing Tobs, indicating an increasing ac- 
cretion rate. So it may be that a particular range in M 
is required to maintain the local Eddington effect. But 
without a clear theoretical answer to these questions, we 
consider a number of intriguing possibilities for the disk 
evolution. 

First, we can take the results of our phase-resolved 
spectroscopy at face value and use our disk parameters to 
infer the disk mass accretion rates for single- and double- 
peaked oscillations (FigurcO . We calculate t he accretion 
rate using Equation 4 of Zimmerman et al.l (2005') . To 
facilitate a visual comparison, we also show the double- 
peaked M vertically normalized to match the mean value 
for the single-peaked cycle (green dashed line). The ac- 
cretion rate is higher on average in the double-peaked 
cycle (M w 7.4 x 10^® g s"^ versus M w 5.9 x 10^^ g 

s~^), but the M wave is wider and has a smaller abso- 
lute amplitude. As an aside, this higher accretion rate 
implies a larger unstable region in the disk, which could 
explain why the cy cle period is slightly lar ger for double- 
peaked cycles (e.g. |N^akshIn£t^n[2003). Considering 
the last seconds before/during the collapse of the disk, 
we see that M is sharply peaked in the single-peaked cy- 
cle, so that for a fixed phase-averaged accretion rate (i.e. 
comparing the normalized accretion rates), more matter 
accretes earlier in double-peaked oscillations. 

But this explanation is only valid to the extent that the 
accretion disk dominates both the mass flux and the X- 
ray emission. If there are outflows from very small radio 
or if there is a second component, e.g. the corona, with 
a non-zero accretion rate or some difference radiative ef- 
ficiency, then Figure [7] is only an approximation. We 
note that our accretion rates are calculated from the disk 
prior to Compton scattering. For this reason, the contri- 
bution of the corona could be further involved if there is 

* Chandra observations of GRS 1915-1-105 indicate the pres- 
ence of wind s, but th ese are typically driven from the oute r disk 
UMiller et all 2008 : Neilsen fc Leell200gl : INeilsen et aI.|[20TTa| [hD. 



significant mass transfer between it and the disk, as sug- 
gested by iJaniuk fc Czernvl (l2005f) a n d verified by our 
fits he re and in INeilsen et al.l ()2011al ). IJaniuk fc Czernvl 
(|2005h consider several different prescriptions for mass 
transfer between the accretion disk and the corona (i.e. 
dependent on the luminosity or the accretion rate in 
the disk). It could be that single- and double-peaked 
oscillations transfer mass differently, and that this con- 
trols the evolution of the local Eddington process in the 
disk. As th e corona has a stabilizing effect on the disk 
(jJaniuk fc C zcrnv 2005) , it may also be that the coronal 
geometry differs between different types of cycles, leading 
to different plunge behavior. 

Figure [7] is also only accurate insofar as the disk can 
be approximated as having a constant radiative effi- 
ciency, optical depth, spectral hardening, and a con- 
stant/uniform accretion rate (in a given phase bin) 
in the inner disk. The accuracy of these approxima- 
tions must be determined from simulations that cal- 
culate and fit spectra at each phase of the cycle and 
compare to the physical b ehavior of the ac c retion disk 
in the simulations. Sin ce iNavakshin et al.l (|2000l ) and 
IJaniuk fc Czernvl (|2005l ) invoke strong variations in the 
disk viscosity and surface density, it may be that our 
approximations here are insufficient for a complete un- 
derstanding of the evolution of the disk. On the other 
hand, it seems clear that the phase dependence of the 
accretion rate (whether as shown in Figure |7| or other- 
wise) is probably the primary factor in determining the 
strength and duration of the X-ray pulses, and may in 
some yet-to-be-determined way influence the duration of 
the local Eddington effect. INavakshin et al.l (|2000[) note 
that the cycle period and profile depend on the mass ac- 
cretion rate in a non-linear way, and it would be very 
interesting to see if future simulations including a local 
Eddington ef fect and mass transfer b etween the disk and 
corona (after lJaniuk fc Czernvl2005l) could reproduce the 
peak structure of the observed cycles. 

The need for future theoretical studies is also high- 
lighted by a recent spectral study o f a long BeypoSAX 
observation of the heartbeat state (jMineo et al.l 120121. 
Dividing the cyc les into 5 characteristic intervals, 
iMineo et al.l (|2012D followed the disk and corona during 
the cycle, and the y observed behavior co nsistent with our 
present work and INeilsen et all (|2011al) : the disk radius 
increases at constant temperature and the coronal elec- 
tron temperature drops when the disk plunges inwards 
and the disk temperature spikes (i.e. in the hard pulse). 
However, their results suggested a decrease in the coro- 
nal optical depth during the hard pulse, which they in- 
terpreted as a condensation of the corona (in contrast to 
the sharp i ncrease in the corona l optical depth reported 
here and in INeilsen et al.ir2011al and interpreted as a re- 
sult of mass ejection processes). Presently, it is unclear 
if this difference is due to their wider energy coverage 
(~ 4x), our superior phase resolution (~ lOx), or dif- 
ferences in our models, i.e. simpl versus compPS (which 
includes non-thermal electrons). But the similarities in 
our results are highly encouraging, and there is hope that 
new theoretical studies will allow us to converge on a 
complete understanding of the p state in the near future. 
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GRS 1915+105 exhibits a fascinating array of variabil- 
ity in its X-ray lightcurve. It is believed that many of 
its variability classes are related to limit cycles of ac- 
cretion, advection, and ejection, but the origin of these 
variability classes is only understood to the degree that 
we understand their diversity. This is particularly true 
for the p state, which is similar to several other variabil- 
ity classes, and which displays striking variations from 
observation to observation, particularly with respect to 
the numb er of peaks per cycle in the X-ray lightcurve 
(e.g. Mas saro et al.ll2010t this work). 

In a previous paper (Neilse n et al.l l2011aD , we per- 
formed a detailed, phase-resolved spectral and timing 
study of GRS 1915+105 in a double-peaked instance of 
the p state. Building on our knowledge of this type 
of oscillation, we have presented here a comprehensive 
study of archival RXTE observations of the p-like cy- 
cles. Broadly, our results indicate that the variation 
in the heartbeat state lightcurves can be explained in 
the context of two main pulses, one soft and one hard, 
whose relative strengths, duration, and separation can 
vary smoothly over time. This explanation appears to 
apply whether each cycle exhibits one pulse, two pulses, 
or some intermediate morphology. Furthermore, phase- 
resolved spectra from the RXTE PCA clearly indicate 
that oscillations with very different profiles exhibit re- 
markably similar spectral evolution. We conclude that 
a common set of dominant accretion processes produces 
the diverse lightcurves. 

In this context, what is particularly remarkable about 
the heartbeat state is the complex interaction between 
the disk, its mass accretion rate, and the radiation it pro- 
duces. The oscillation is a process set in motion by a high 
external accretion rate, in which an unstable disk pcriod- 
ically sends wave s of mass inwards (e.g. Taam et al. 1997; 
iNavakshin eFall [2OOOI : iJaniuk fc Czmwl [2005t ). Re- 
sponding to the increasing influx of matter, the disk pro- 
duces enough radiation to slowly expel some fraction of 
the inflow, driving the inner edge of the disk away from 
the black hole. But eventually, the disk's ability to com- 
pensate is overwhelmed; when the accretion rate remains 
high, the excess matter in the inner disk plunges inwards 
as the temperature and viscosity spike. But even in these 
final moments of the cycle, radiation pressure still man- 
ages to eject some mass from the disk. These are the 
common, interconnected processes that lead to the p-like 
cycles in GRS 1915+105. 

But there is another process, the accretion disk wind, 
whose ubiquity has yet to be confirmed (see our up- 
coming Chandra, RXTE, Gemin i, and EVLA study ; 
Neilsen et al., in preparation). In lNeilsen et al.l (|2011aD 
we showed that the bright X-ray pulses from the in- 



flow drive a massive, ionized wind off the outer disk 
[R < 10^^ cm). This wind can drain as much as 95% 
of the external mass supply from the accretion flow, 
suggesting that the companion star is actually supply- 
ing matter at a rate well in excess of Eddington (since 
matter reaches the black hole at roughly the Edding- 
ton rate). Mass ive winds have a l so been seen in o ther 
black holes (e.g. iKing et al.l 120121: iPonti et al.ll2012l and 
references therein ) and neutron stars (e.g. GX 13+1; 
lUeda et al.l 120041 ). Eventually, this massive wind de- 
pletes the outer disk, creating a mass deficit. The heart- 
beat oscillations continue unabated until, days or weeks 
later, this deficit propagates inwards dShields et al.|[l986l : 
iLuketic et al.ll201u iNeilsen et al.ll2011aD . simultaneouslv 
suppressing the thermal- viscous instability, the local Ed- 
dington effect, radiation-pressure-driven mass ejection, 
and the accretion disk wind itself. 

We end on an exciting note, for it is finally clear 
that GRS 1915+105 is not alon e. In a new out- 
burst detected by Swift/BAT (jKrimm et al.l 120111) 
and followed up with frequent RXTE PCA obser- 
vations, the black hole candidate IGR J17091-3624 
was discovered to exhibit heartbeats and several other 
of GRS 1915+ 105's variability classes in its X-ray 
lightcurve (e.g. iPahari. Yadav. Sz Bhattacharvval 120111 : 
lAltamirano et al.ll2011[ ). But if GRS 1915+lQ5's behav- 
ior is made possible by a high mass supply rate from its 
subgi ant companion (Greiner, Cu bv. fc McCaughreanI 
120011 ). is the presence of two active "GRS 1915+105"s in 
our Galaxy consistent with models of stellar and binary 
evolution? Thus, although at a much lower X-ray flux, 
future studies of IGR 17091-3624 will place strong con- 
straints on the origin of accretion instabilities like those 
in GRS 1915+105. 
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